Fluctuation theorems for capacitively coupled electronic currents 
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The counting statistics of electron transport is theoretically studied in a system with two capaci- 
tively coupled parallel transport channels. Each channel is composed of a quantum dot connected 
by tunneling to two reservoirs. The nonequilibrium steady state of the system is controlled by 
two affinities or thermodynamic forces, each one determined by the two reservoirs of each channel. 
The status of a single-current fluctuation theorem is investigated starting from the fundamental 
two-current fluctuation theorem, which is a consequence of microreversibility. We show that the 
single-current fluctuation theorem holds in the limit of a large Coulomb repulsion between the two 
parallel quantum dots, as well as in the limit of a large current ratio between the parallel channels. 
In this latter limit, the symmetry relation of the single-current fluctuation theorem is satisfied with 
respect to an effective affinity that is much lower than the affinity determined by the reservoirs. 
This back-action effect is quantitatively characterized. 



I. INTRODUCTION 



Away from equilibrium, fluctuating currents flow 
across small open quantum systems such as quantum 
dots exchanging electrons with reservoirs. Advances in 
nonequilibrium statistical mechanics have shown that 
the current fluctuations obey symmetry relations fol- 
lowing from microreversibility and known as fluctuation 
theorems^— They have been proved in different contexts 
and, especially, for open quantum systems and the full 
counting statistics of electron transport In this con- 
text, fluctuation theorems relate the probabilities of op- 
posite random values of the currents to the potential dif- 
ferences driving the mean values of the currents. In elec- 
tronic circuits, these potential differences play the role 
of thermodynamic forces also called affinities; 16 ! 17 Fluc- 
tuation theorems hold in nonlinear transport regimes, in 
particular, for the description of the Coulomb drag effect 
in capacitively coupled quantum dotsJ^ 

Remarkably, modern technology is able to perform 
the bidirectional counting of single-electron transfers in 
quantum-dot circuits, allowing the experimental verifica- 
tion of the fluctuation theorem* 1 ^ In these experiments, 
the quantum-dot (QD) circuit is monitored by a parallel 
circuit made of a quantum point contact (QPC). Because 
of electrostatic interactions, the electronic occupancy of 
the quantum dots modifies the mean value of the QPC 
current, enabling the measurement of the QD electronic 
state in real time. The surprise has been that, within 
experimental error, the bidirectional counting of the QD 
current obeys the symmetry relation predicted by the 
fluctuation theorem but with respect to an affinity about 
one order of magnitude smaller than the potential dif- 
ference driving the QD circuit* 1 ^ This discrepancy has 
revealed the importance of the interaction between the 
QD and QPC circuits. Indeed, the QPC current is typ- 
ically 10 7 -10 8 times larger than the QD current in such 
experiments so that the QPC can act as a quasi-classical 
detector measuring the quantum state of the QDs. As a 
consequence, the whole system composed of the two par- 



allel circuits is quite far from equilibrium and the shot 
noise in the QPC current has a significant back action 
onto the small QD current. In Ref. Q11 this back action 
was analyzed in terms of the so-called P{E) theory^ by 
fitting experimental data to a simple Lorentzian in order 
to take into account the global effect of the QPC noise 
onto the QD tunneling rates. In Ref. [2l|, a stochastic 
model was proposed with extra QPC states besides the 
QD states. Both approaches leave open the fundamen- 
tal understanding of the back action in terms of the mi- 
croscopic Hamiltonian describing the interaction between 
the QD and QPC circuits. 
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FIG. 1: Schematic representation of two quantum dots in 
parallel. Each quantum dot is coupled to two reservoirs of 
electrons. Moreover, both quantum dots influence each other 
by the Coulomb electrostatic interaction. 

In the present paper, we address this issue by con- 
sidering a system composed of two capacitively coupled 
parallel transport channels, each containing a single QD 
in contact with two electron reservoirs! 14 ' 22 ! 23 The two 
parallel transport channels are only coupled by the elec- 
trostatic Coulomb repulsion between the electrons occu- 
pying the two QDs so that there is no electron trans- 
fer between both channels (see Fig. [T]). The two cur- 
rents flowing in parallel across this device are driven by 
two affinities defined by the potential differences on both 
QDs. In the following, the circuit No. 1 plays the role of 
the QD and the circuit No. 2 the role of the QPC. The 
fluctuations of the two currents obey a bivariate fluctu- 
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ation theorem, which is the fundamental consequence of 
microreversibility. Under general conditions, this two- 
current fluctuation theorem does not imply the existence 
of a single-current fluctuation theorem for the QD cur- 
rent monitored by the secondary QPC circuit. Here, 
we show that the single-current fluctuation theorem only 
holds in the limit where the current in the secondary cir- 
cuit is much larger than the one in the main circuit (or 
vice versa). However, the symmetry of the single-current 
fluctuation theorem does not hold with respect to the 
potential difference on the main circuit but to an effec- 
tive affinity which strongly depends on the electrostatic 
interaction between both circuits. In this way, our anal- 
ysis provides an understanding of these features in terms 
of the basic parameters of the system Hamiltonian and 
clearly shows that the effective affinity can vary by one 
order of magnitude or more due to the back action of the 
QPC onto the QD. 

Furthermore, our analysis leads to the evaluation of the 
entropy production in the electronic device. The fluctu- 
ation theorem has for consequence the non-negativity of 
the entropy production and is thus compatible with the 
second law of thermodynamics. The directionality due to 
the nonequilibrium driving of the device is characterized 
by the probability distributions of the current fluctua- 
tions, by the mean values of the currents, and also by 
the entropy production. The analysis based on the fluc- 
tuation theorem allows us to understand the connections 
between these complementary and fundamental aspects 
of such nonequilibrium electronic devices. 

The paper is organized as follows. In Sec. [Til the 
Hamiltonian model is presented and the master equation 
ruling the occupancies of the QDs is derived for QDs 
weakly coupled to the reservoirs within the Markovian 
and secular approximations. Section |Hl] is devoted to 
the full counting statistics of the two interacting currents, 
for which the fundamental fluctuation theorem is estab- 
lished. Moreover, the connection between the fluctua- 
tion theorem and the entropy production of the device 
is discussed. In Sec. IIV1 we first consider the limit of 
a large Coulomb repulsion between the QDs, in which 
case a single-current fluctuation theorem is obtained but 
without modification of the effective affinity contrary to 
the experimental observation. In Sec. [Vj the limit is then 
considered where the current in one circuit is much larger 
than the one in the other circuit. It is in this limit that 
the single-current fluctuation theorem is obtained with 
an important modification of the effective affinity with 
respect to which the symmetry relation of the single- 
current fluctuation theorem holds. In Sec. IVI1 these 
effects arc numerically demonstrated with the model for 
parameter values corresponding to typical experimental 
conditions. We analyze the dependence of the effective 
affinity on the parameters of the Hamiltonian model and, 
especially, on the electrostatic interaction between both 
circuits. Conclusions are drawn in Sec. IVIII 



II. CAPACITIVELY COUPLED PARALLEL 
TRANSPORT CHANNELS 

A. The Hamiltonian 

The vehicle of our study is the Hamiltonian model con- 
sidered in Ref. 22;. Each transport channel (a = 1 or 
a = 2) is composed of one quantum dot with a single 
energy level e a for the electron. This level is either occu- 
pied or empty and the spin degree of freedom is ignored. 
Moreover, the quantum dots are capacitively coupled by 
electrostatic repulsion if both are occupied. This electro- 
static repulsion is taken into account by an Anderson- 
type term with the parameter U. The parameter U is 
thus the energy contribution of the Coulomb repulsion 
when both quantum dots are occupied by an electron. 
The system Hamiltonian is therefore given by 

H s = e 1 d\d 1 +e 2 d\d 2 + Ud\d l d\d 2 (I) 

where d a and dJ a denote the annihilation and creation 
operators of an electron on the QD labeled by a = 1,2. 
This Hamiltonian is diagonalized in the four-state ba- 
sis {|00}, 1 10), 1 01} , |11)} with the corresponding energy 
eigenvalues {0, ei, e 2 , t\ + e 2 + U}. 

Each QD is in tunneling contact with two reservoirs 
on its left- and right-hand sides (see Fig. []}. The system 
has thus four reservoirs j = 1L, 1R, 2L, 2R, which are 
denoted as j = ai by the label a = 1, 2 of the channel 
and the side i = L,R where the reservoir stands. The 
Hamiltonian of all the reservoirs can be expressed as 

H R = J2 H i ( 2 ) 

j 

in terms of the Hamiltonians of the individual reservoirs, 
which are defined as 

H i = Yl e 3 k c )k c 3k ( 3 ) 

k 

where Cjk and cj fc denote the annihilation and creation 
operators of electrons in the corresponding states. The 
reservoirs are supposed to be much larger than the sys- 
tem itself so that the eigenvalues {tjk} of each reservoir 
form a very dense spectrum which is quasi continuous and 
characterized by a density of states Dj(e) = J2 k S(c—£jk)- 
The operator giving the electron number in the reservoir 
j is furthermore defined as Nj — J2k c ]k c jk- 

The tunneling Hamiltonian establishing the interaction 
between the QDs and the reservoirs has the form: 

HSR = ^ ^ X! taik d a°aik + H. C. (4) 
Q =l,2i=L,R k 

where we have here specified the channels and the reser- 
voirs by writing j = ai. The effect of the electrostatic 
interaction on the energy barriers between the quantum 
dots and the reservoirs could be taken into account by 
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including corresponding capacitances, as considered in 
Ref. This effect is neglected in the Hamiltonian model 
of Ref. Hil that we here consider. 

Finally, the total Hamiltonian is defined as the sum: 



H — Hs + i?R + Hsr 



(5) 



We notice that the electron number operators of each 
transport channel 

N a = dUa + ]T ]T c f mk c mk a = 1, 2 (6) 

i=L,R k 

separately commutes with the total Hamiltonian 

[H, Ni] = [H, N 2 ] = (7) 

so that the electron number is conserved on each trans- 
port channel and there is no electron exchange between 
the channels. In contrast, the number operators of the 
reservoirs Nj = N ai with a = 1, 2 and i = L, R do not 
commute with the total Hamiltonian unless the tunneling 
amplitudes are equal to zero. 



B. The master equation 

Initially, the reservoirs are in grand-canonical equilib- 
rium states characterized by the chemical potentials fij 
with j G {1L, 1R, 2L, 2R} and a uniform temperature 
T. We denote by j3 — (feeT) -1 the inverse temperature 
with the Boltzmann constant k^. On the other hand, the 
QDs are in an arbitrary statistical mixture /?s(0). More- 
over, a measurement is performed at the initial time that 
determines the numbers mi and 1112 of electrons in the 
reservoirs j — 1L and j — 2L. Consequently, the ini- 
tial density matrix of the total system has the factorized 
form: 



Pi 



iV, 



'iV 1L ,mi 0jv 2Lj m2 



(8) 

where Hj denotes the partition function of the grand- 
canonical ensemble for the reservoir j. The Kronecker 
symbols £jv,m take the unit value if N = m and zero oth- 
erwise and they thus play the role of projection operators 
on states with a fixed number of particles. Thereafter, 
the density matrix of the total system evolves in time 
according to the Landau- von Neumann equation 



idtp mi m 2 {t) = [H, p mim2 (t)] 



(9) 



in units where h = 1. The following normalization con- 
dition is satisfied by the initial density matrix ([5]) and 
preserved by the time evolution ([9]): 



^2 Tr Pm im2 (t) = 1 



(10) 



where Tr denotes the trace over all the degrees of freedom 
of the total system. 



Since we are interested in the occupancies of the QDs 
and the numbers of electrons transferred between the 
reservoirs, we focus on the probabilities p VlV2 {n\, n-i) that 
the QDs are in the quantum states {|^i^2}} with the oc- 
cupancies v\ = 0, 1 and v 2 = 0, 1, while n\ electrons have 
been transferred from the reservoir j = 1L to the first 
QD and n 2 electrons from the reservoir j = 2L to the 
second QD between the initial time t = and the time t. 
These probabilities can be defined in terms of the density 
matrix of the total system according to 

(t) \ V 1 V 2 ){V 1 V2\ 

mi ,7712 

x ^JViL,mi— m &N 2 - L ,m 2 ~n 2 (H) 

which results from a second measurement at time t count- 
ing the numbers of transferred electrons*^ 

We suppose that the two quantum dots are weakly 
coupled to the reservoirs by small enough tunneling am- 
plitudes {tjk} so that we may carry out the Born per- 
turbative approximation on the Landau-von Neumann 
equation up to second order in the tunneling amplitudes. 
We use the secular (or rotating wave) approximation and 
we take the Markovian approximation . 11 ' 24 " — As a con- 
sequence, the charging and discharging transition rates 
of the QDs are respectively given by 



a j — Ljfj 

in terms of the Fermi-Dirac distributions 

1 



f= — 

1 

1 _|_ ^{ej+U-^j) 



fj 



(12) 
(13) 
(14) 
(15) 



(16) 
(17) 



where ej = e a for j = ai. The rate constants are given 
at the second order of perturbation theory by 



k 

= 2TT\t 3 {e j )\ 2 D J {e 3 ) 
f 3 = 2TT^\t jk \ 2 5(e 3 + U - t jk ) 



(18) 



= 2^( £j + U) \ 2 D {e J+ U) (19) 

where the quantities tj(e) are the tunneling amplitudes 
as a function of energy and Dj (e) the density of states of 
the reservoir j. 

The total system is characterized by two sets of time 
scales: 

(1) The correlation times of the reservoirs: The cor- 
relation time of the reservoir j can be estimated as 
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Ae ■ in terms of the width Ae, of the function 



3 3 

giving the charging rate a,j(e) — 2ir\tj(e)\' 2 Dj(e)fj(e) ver- 
sus the energy e. 

(2) The relaxation times induced by the electron ex- 

(R) — l 

changes with the reservoirs: rj ~ Tj . 

In consistency with the assumption of weak coupling, 
we suppose that the correlation times are much shorter 
than the relaxation times and that the secular approxi- 
mation is performed by averaging the equation of motion 
over a time scale At which is intermediate between both 



which justifies the use of the Markovian 
approximation . 24 ' 25 Moreover, the Lamb shifts of 
the QD energy levels are neglected for simplicity. 

Accordingly, the master equation for the probabilities 



p(ni,n 2 ) = 



/ poo(m,n 2 ) \ 
Pio(m,n 2 ) 
Poi(ni,n 2 ) 

\ pn(ni,n 2 ) I 



(21) 



rj C) « At « rf > 



(20) 



takes the form 



d t p(ni,n 2 ) = (Li + L 2 ) ■ p(n 1 ,n 2 ) 



with the matricial operators: 
L 1 = 

and 



(22) 



/ -aiL-aiR b lh E^ + b 1R \ 

aiLE^+am -6il - &m 

-aiL - Sir biL_E+ +J>m 

\ aiL El + ai R -6 iL - &m / 



/ -a 2L - a 2R b 2L E+ + b 2R ^ ^ 

-a 2 L - a 2 R b 2 h E^ + b 2 R 



(23) 



V 



«2L ^2 + °2R 






02L ^2" + "2R 



— b 2 L — b 2 R, 








(24) 



32L 



&2R / 



where the step operators 



E d , 



exp ± 



dn r 



(25) 



III. THE TWO-CURRENT FLUCTUATION 
THEOREM AND ITS CONSEQUENCES 



The cumulant generating function and the 
affinities 



increase or decrease the numbers n a of transferred elec- 
trons 

Etf{n a ) = f(n a ± 1) (26) 

when applied on any function /(n a )i2I 

We notice that the occupancy probabilities irrespective 
of the numbers of transferred electrons defined as 

-foo +00 

p uiu 2 = ^2 ^2 Pu 1 u 2 (ni ) n 2 ) (27) 

m= — 00 n-2— — 00 

obey the master equation obtained by replacing the step 
operators (|25p by unity, E^ = 1, in the matricial opera- 
tors (HI) and 423} . 



In order to perform the counting statistics of the elec- 
trons transferred from the left reservoirs to the quantum 
dots, we introduce the cumulant generating function of 
the currents in terms of the counting parameter X a of the 
corresponding transport channel: 

Q(Ai,A 2 )= lim -- In (exp(-Aini - \ 2 n 2 )) t (28) 

t— >oo t 

where the average 

(X)= V^u 2 {n x ,n 2 )X (29) 

is taken with respect to the probability distribution, 
which is the solution of the master equation (f22j) at the 
time t. 

We notice that the cumulant generating function (|2"5)l 
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is given as the leading eigenvalue of the following eigen- 
value problem: 

L-v = -Qv (30) 

where 

L = c- A " + L 2 ) c+ An = Li + L 2 (31) 

is a four-by-four matrix with real elements that de- 
pend on the counting parameters A. Since the functions 
exp(A • n) are the eigenfunctions of the step operators 



(1251) . the matrices l_i and L 2 are obtained by the follow- 
ing substitutions 

Et -> e ±A ° (32) 

in Eqs. (1231) and (|2~4")) , as can be checked by a straight- 
forward calculation. 

The four- by-four matrix L = L(A) obeys the symmetry 

M 1 L(A) M = L(A A) T (33) 

with 



M 



/l 


V o o 






(C2— A»2r) 









(ei+e2 + !7-AiiR— M2r) 



(34) 



and the affinities A — (A\, A 2 ) defined by 



or 



Ai = In = in 



A 2 = In 



hham 

Q2L&2R. 
b 2 LO, 2 R 



In 



^2L«2R 



/3 (Mil -Mir) (35) 
= /3 (M2L - M2r) (36) 



We notice that these affinities can also be obtained by 
using Schnakenberg graph analysisi&2£ These quantities 
are the two independent thermodynamic forces able to 
drive the system away from equilibrium. The fact that 
there exists only two independent affinities although the 
system contains four reservoirs is due to the existence 
of the two constants of motion (0 given by the particle 
numbers in the two transport channels. 

If the system was fully connected, only the total par- 
ticle number would be a constant of motion and there 
would exist three independent affinities. More generally, 
a system composed of r reservoirs and partitioned into 
c disconnected but interacting transport channels has c 
constant particle numbers and can be driven away from 
equilibrium by r — c independent affinities. Here, r = 4 
and c — 2 so that there is only r — c = 2 independent 
affinities. 

As aforementioned, the cumulant generating function 
is given by the leading eigenvalue of Eq. (|3"0")) , i.e., by the 
smallest root of the quartic characteristic polynomial: 



det(L 







(37) 



of the four- by- four matrix pip . Therefore, the symme- 
try (j3"3"|) implies that the cumulant generating function 
also obeys this symmetryj&£ In this way, the fundamen- 
tal result is proved that the cumulant generating function 
satisfies the fluctuation theorem: 



Q(X) = Q(A - A) 



(38) 



Q(\ u \ 2 )=Q{A l -\ l ,A 2 -\ 2 ) 



(39) 



in terms of the affinities A = [A\ , A 2 ) given by Eqs. (|3"5"|) 
and ([56). 

An alternative expression of this fluctuation theorem 
is that the probability 



p{m,n 2 )= y p VlV2 {m,n 2 ) 



vi,v 2 



(40) 



for the transfer of n\ particles in the circuit No. 1 and 
n 2 particles in the circuit No. 2 during the time interval 
t obeys 



p{n u n 2 ) 
p(-ni, -n 2 ) 



exp(Aini + A 2 n 2 ) for t — > +oo 



(41) 

Indeed, this expression implies Eq. (131?)) using the defi- 
nition (I28p of the cumulant generating function with the 
average (j2U|) £ili 

In general, this two-current fluctuation theorem does 
not imply any single-current fluctuation theorem unless 
specific conditions are satisfied either by construction, 29 
or in some particular limit, which is here the case as 
shown in the following sections. 



B. The average currents and the response 
coefficients 



The average values of the particle currents are given in 
terms of the generating function according to 



J cy — 



dQ 

8X a 



(42) 



A=0 
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for a = 1,2. As shown in Appendix \K\ these currents 
can be expressed in terms of the probabilities (|2T|) of 
the four QD states in the nonequilibrium steady state 
corresponding to the affinities (|3"5)) and (155)) according 
to: 



Jl — OlL-PoO 
J 2 = 02L-Po0 



h h P W 
&2L-Poi 



- CtlL-Poi 

a2L-Pio 



&1L-P11 

62LP1I 



(43) 
(44) 



These currents are nonlinear functions of the affinities, 
which can be expanded in powers of the affinities in order 
to identify the linear and nonlinear response coefficients: 

J a = J a (A 1 ,A 2 ) 

= ^L a ^A & + -^Ma^ApAy 



0,1 



(45) 



As a consequence of the fluctuation theorem (|39p . the 
linear response coefficients L a ^ are given in terms of the 
second derivatives of the generating function with respect 
to the counting parameters and they thus obey the On- 
sager reciprocity relations: 



L a ,p — Lft_ n — — — 



1 d 2 



2 dX a d\r- 



A=0.A=0 



(46) 



Similar relationships have been established for the non- 
linear response coefficients.— 

The average currents as well as the linear response co- 
efficients can be calculated in terms of the characteris- 
tic determinant (|37p of the matrix (|31[) as shown in Ap- 
pendix [A] By using Eq. (|A10I) , the Onsager coefficient 
turns out to be proportional to 

Liz oc (r 1L f ir - f uPih,) (r 2L f 2R - f 2L r 2 R) (47) 

In general, the Onsager coefficient is thus non vanishing 
and there is a phenomenon of Coulomb drag according 
to which a current may be induced in a circuit at equilib- 
rium if the other circuit is out of equilibrium, as shown 
in Ref . El 

However, the Onsager coefficient vanishes under the 
condition that the rate constants of one circuit do not 
depend on the Coulomb repulsion parameter U . In this 
particular case, there is no Coulomb drag because 



Ji(0,A 2 ) = 



and 



J a (4t,0) = Q 



if r, =1^ 



(48) 

This property is also proved in the Appendix [A] Equa- 
tion (|4"5|) implies the vanishing of the Onsager coefficient 
as well as the nonlinear response coefficients allowing the 
coupling of one current to the affinity of the other circuit: 

ii.2 = M\,22 = ^1,222 = • • • = and 

£2,1 = M 2 ,u = N 2tlu - • • • = if Tj = f d (49) 

Nevertheless, these coefficients do not vanish in general. 



C. The entropy production and the energy 
dissipation 

A further consequence of the fluctuation theorem 
is that the average currents (|4"2"|) obeys the second law of 
thermodynamics according to which the entropy produc- 
tion is always non-negative: 



- — ±- = A 1 J 1 + A 2 J 2 >0 
kr at 



(50) 



where fee is Boltzmann's constanti&ii 

The power dissipated in each circuit is defined as the 
product of the voltage V a by the electric current I a — eJ a 
where e is the electric charge of the particle: ±l a = V a I a , 
with a — 1,2. Since the affinities are related to the 
voltages by 



A n 



(51) 



we have that the dissipated power in the circuit a is given 
by 



II a — /sb T A a J a 



(52) 



and the entropy production is thus proportional to the 
total dissipated power: 



djS 



(III +n 2 ) > 



(53) 



Therefore, the entropy production of the system charac- 
terizes the energy dissipation of the quantum measure- 
ment performed on one QD by the current flowing in the 
other circuit p laying the role of the QPC in the experi- 
ments of Rcfs. 18.31. We shall evaluate these quantities 
under such specific conditions in the following sections. 



IV. THE LARGE COULOMB REPULSION 
LIMIT 

In the limit where the Coulomb repulsion between both 
QDs is large, the coupling parameter U takes large values 
so that the charging rates of a second electron on the two 
QDs vanish: 

tij = for j = 1L, 1R, 2L, 2R if U = oo (54) 

As a consequence, the probabilities (|TT|) and (|2T[) that 
the system is in the fourth state |11) also vanish: 



Pn(ni,n 2 ) 







and 



Pu 







if U = oo 

In this limit, the occupancy of one QD is stopping the 
current in the other QD. For instance, the average cur- 
rent in the secondary circuit (|44|) has two contributions 
depending on the occupancy of the QD No. 1: 



■h 



1^=1 



(56) 



7 



However, the contribution when the QD No. 1 is occupied 
is vanishing 



=1 = a 2L Pio - &2LP11 =0 if U = 



(57) 



since a 2 L = according to Eq. ((54]) and P n = because 
of Eq. (|55|) . Therefore, the secondary circuit has a non- 
vanishing current only when the QD No. 1 is empty and 
vice versa. 

The cumulant generating function can be obtained by 
considering the three-by-three matrix obtained by remov- 
ing the fourth row and column from the matrix (|31[) . In 
this case, the characteristic determinant (|37p depends on 
the counting parameters only in the following combina- 
tions: 



am,b 1L e Xl 
a 2 R&2Le A2 



(58) 
(59) 



which remain invariant under the independent substitu- 
tions Ai — > A\ — \\ and/or A2 — > A 2 — A2 with the affini- 
ties (1331) and (j3"6")l . Consequently, we obtain the symme- 
try relations: 

Q(Ai,A 2 ) = Q(A 1 -Ai,A 2 ) 

= <9(Ai,A 2 -A 2 ) = Q(^i-Ai,A 2 -A 2 ) (60) 

if U = 00, which implies the single-current fluctuation 
theorem: 



Q(Ai,0) = Q(Ai-Ai,0) if U 



00 



(61) 



but with the unmodified affinity (|35p . Therefore, this 
limit cannot explain the modification of the affinity ob- 
served in the experiments reported in Ref. Il8l . 



V. THE LARGE CURRENT RATIO LIMIT 

In counting statistics experiments! 18 ' 31 , the QPC 
which is used to observe the occupancy of the QD carries 
a current which is typically much larger than the current 
in the QD by a huge factor 10 7 -10 8 . If the QPC is taken 
as the circuit No. 2 in the present model, the rate con- 
stants of that circuit are much larger than the ones of 
the circuit No. 1: 



Til, Tir, Til, Tir, <C r 2 L, r 2 R, r 2 L, T 2 r 



(62) 



Under such circumstances, the relaxation times ~ 
r^ 1 of the circuit No. 1 are much longer than the relax- 



ation times T 2i ~ of the circuit No. 2 and the mon- 
itoring of the slow circuit by the fast one is performed 
over a time scale At such that 



r|f > « At « rj 



(R) 



(63) 



instead of the time scale (|2U|) . 

Our aim is here to obtain the cumulant generating 
function for the counting statistics in the sole circuit No. 1 



without measuring the current in the fast circuit No. 2, as 
it is the case in Refs. Il8ll3~ll . This amounts to consider the 
two-current generating function (f2"5| for A2 = 0. Accord- 
ingly, we focus on the time evolution of the probabilities 
defined by 



1^2—0,1 ri2 — — o 



Pv 1 y 2 {n\,n 2 ) 



(64) 



Since the electron transfers in the circuit No. 2 are 
much faster than in the circuit No. 1, the circuit No. 2 
can be supposed to be in a stationary state during the 
whole period when the circuit No. 1 is in a given state. 
Such stationary states conditional to the state v\ of the 
QD No. 1 are obtained by finding the zero eigenvectors of 
the transition matrix (l24l) with — 1 • The conditional 
probabilities P v% \ Vl that the QD No. 2 has the occupancy 
v 2 provided that the QD No. 1 is in the state V\ are given 
by 



M)|0 — 

A|o = 

P \l = 



a 2 + 


l>2 


a 2 




a 2 + 


t>2 


h 




a 2 + 


h 


0,2 




a 2 + 


h 



with 



"2 
b 2 
a-2 



- 02L + »2R 
&2L + &2R 
: a 2 L + 02R 



bo = b 



2L 



J2R 



(65) 
(66) 

(67) 
(68) 



(69) 
(70) 
(71) 
(72) 



Under the conditions (l6"2l , the probability that the sys- 
tem is in the state \viv 2 ) and that ri\ electrons have been 
transferred in the circuit No. 1 factorizes into the proba- 
bility (|6"4"]) and the probability of the occupancy v 2 of the 
QD No. 2 conditioned to the occupancy v\\ 



(73) 



Substituting these relations into the master equation 
(|22[) and summing over n 2 and v 2 , we get the master 
equations for the probabilities p Vl (n{) as follows: 

dtPo(ni) = -fa L + an) Po(»i) + (h,E^ + b R jpi(n 1 ) 

(74) 

dtpi(ni) = (a L E± + aj^j po(m) - (b^ + &R.J Pi(ni) 

(75) 



where 



a-L = aiL-Po|o + aiL-Pi|o 
or = aiR-Fb|o + am-Pi|o 
bh = &il-Po|i + bihPi\i 
bR = &1R-P0I1 + hnPm 



(76) 
(77) 
(78) 
(79) 
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are the charging and discharging rates of the first quan- 
tum dot averaged over the conditional stationary prob- 
abilities of the second quantum dot. The master equa- 
tions (T74")) - (rT5l rule the process in the slow circuit No. 1 
as monitored by the fast circuit No. 2 over the time scale 

(EH). 

Taking a solution of the form p Vl (ni) ~ exp(Aini — 
Qt) for Eqs. f71 |) -([75 ]) . the cumulant generating function 



J 



(p8|) with A 2 = has thus for approximation the leading 
eigenvalue of the matrix 



6 L e +Al + & E 

QLe'^+OR "L -"II 



1 _ 1 -ah - a R o L e 1 — -f o R 
L - 1 - -br, - h 



which is given by 



Q(Ai,0) 



«L 



+ a R + b L + b R - \J (a L + a R - b L - b R f + 4 (a L c- A i + a R ) (6 L e+ A i + b R ) 



(80) 



(81) 



in the limit (|62[) where the current in the second quantum 
dot is much larger than in the first one. In this limit, the 
generating function (|81[) obeys the single- current fluctu- 
ation theorem: 



Q(Ai,0) = Q(i 1 -Ai,0) 



(82) 



with the effective affinity for the first quantum dot ob- 
tained as 



A x = In 



or&l 



(83) 



in terms of the averaged rates (1751) - ((75)) . This constitutes 
the main result of the present paper. 

We notice that similar results hold in the other limit 
where the circuit No. 1 is much faster than the circuit 
No. 2 because both circuits have the same structure and 
are symmetrically coupled together through the Coulomb 
repulsion of parameter U in Eq. (TTJ). 

The result ([82]) shows that the generating function of 
the counting statistics in the slow QD No. 1 has the sym- 
metry of a single-current fluctuation theorem under the 
experimental conditions (|62[) but with respect to the ef- 
fective affinity (|53")) . This latter may differ by orders of 
magnitude with respect to the affinity (|35l) driving the 
circuit out of equilibrium. The reason for this modifica- 
tion is the back action of the other circuit to which the 
QD is capacitively coupled. Indeed, the charging and 
discharging rates of the QD No. 1 are averaged over the 
two possible states of the QD No. 2 according to Eqs. 
(fT6")) - (|79[) so that their effective values are modified by 
the back action of the circuit No. 2. This modification 
of the transition rates is reminiscent of the influence of 
environmental noises as described by the P(E) theory.— 

In the following section, the importance of the back ac- 
tion is numerically demonstrated by considering specific 
conditions. 



VI. NUMERICAL RESULTS 

A. Parameter values 

In typical counting statistics experiment^"- 19 the 
affinities take quite large values because the voltages are 
large with respect to the temperature. For instance, the 
voltages Vqd = 300 fxV, Vqpc = 800 fxV, and the elec- 
tronic temperature T = 130 mK are reported in Ref. [l8l . 
In our analogy, we choose the affinities as follows: 



.1, A QD ^ 



A 2 = A, 



QPC 



k B T 
eVQPC 
k R T 



= 70 



(84) 
(85) 



Since the QPC current is reduced by about 10% if the QD 
is occupied, the parameter U of the Coulomb repulsion 
between both QDs can be taken as 



0U = 32.8 



(86) 



Moreover, the QPC current is about 10 7 -10 8 larger than 
the QD current. 

As aforementioned, the role of the QPC is played by 
the circuit No. 2 and the QD by the one of the circuit 
No. 1 in our model. The energy level of the second 
quantum dot is supposed to be in the middle between the 
reservoirs electrochemical potentials and the couplings to 
the reservoirs are chosen symmetric and independent of 
the energy. Under such assumptions, possible parameter 
values are given by 



(87) 





= 25 








PfJ-m 


= 








r iL = 


Tir 


= flL 


= r 1R = 


1 




= 70 








/5M2R 


= 








r 2L = 


r 2 R 


= f 2L 


= r 2 R = 


10 


/?e 2 = 


35 









(89) 
(90) 
(91) 
(92) 
(93) 
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while the level of the QD No. 1 has the energy e±, which 
may take different values in the following numerical cal- 
culations. The correlation times of the reservoirs are sup- 
posed to be short enough for the conditions ([20)) to be 
satisfied in consistency with the perturbative approxima- 
tion. 



B. Stochastic simulations 

The random time evolution of the system can be gen- 
erated by simulating the stochastic jump process of the 
master equation (|22|) with Gillespie's algorithmi 32 i 33 Four 
possible transitions may occur from each of the four 
states. The transition rates are given by Eqs. (fT2"]) - (fT5)) 
with the Fermi-Dirac distributions (| 1 6 1) - (p~7)) and the rate 
constants ([89" ]> -(|92 ]) . 



5 10 7 

I 

4.5 10 7 





1 1 1 




1 i i i i 




I 







10 20 30 40 

time 



FIG. 2: Simulation with Gillespie's algorithm of the QPC 
current in circuit No. 2 measuring the QD occupancy. The 
parameter values are given by Eqs. (|86|l - (|93| l and /3ei = 0. 
The effective affinity of the circuit No. 1 is Ai = 1.17. The 
mean value of the QD current is Ji ~ 0.17 electrons per unit 
time. The mean value of the QPC current is J2 ~ 4.8 x 10 7 
electrons per unit time. The QD is empty (resp. occupied) 
when the QPC current takes the value 5x 10 r (resp. 4.5 x 10 7 ). 

Figure [5] depicts the current in the circuit No. 2 av- 
eraged over a time interval At = 0.01, which is shorter 
than the typical dwell time of the QD No. 1, as required 
by Eq. ([63]) . We see that the current is reduced by about 
10% when the QD No. 1 is occupied, which is in agree- 
ment with the choice for the parameter (|5rJ|. The ratio 
between the mean values of the currents is here given by 
J2/J1 = 2.8 x 10 s , while the ratio of the dissipated pow- 
ers takes the value n 2 /IIi = (A 2 J 2 )/{A 1 J 1 ) = 7.9 x 10 8 . 
Such very large ratios are required in order for the sec- 
ondary current to distinguish between the two states of 
the QD in the primary circuit. Simulations show that the 
fluctuations of the secondary current would be larger for 
smaller values of the current ratio. Thanks to the large 
ratio, the instantaneous occupancy in the circuit No. 1 
can be monitored by the current in the circuit No. 2 over 
the time scale (|rj3")l . which is longer than the time scale 
of the fast circuit No. 2 but shorter than the one of the 
circuit No. 1. 



C. The cumulant generating function and its 
properties 

The cumulant generating function of the current in the 
circuit No. 1 is calculated by the leading root of the char- 
acteristic polynomial ([37)1 of the four-by-four matrix (j3"Tj) 
with A 2 = 0. 

The lack of symmetry of the single-current generat- 
ing function Q(Ai,0) is manifest if the rate constants of 
both circuits are of the same order of magnitude. The 
generating function and its symmetric with respect to 
the effective affinity is depicted in Fig. [3]for r 2s /ri s = 1 
(with s — L, R) and j3U — 30. Here, the effective affinity 
is taken as the non-trivial root of the generating func- 
tion such that Q(Ai,0) = 0. We clearly see that the 
generating function is not symmetric with respect to the 
effective affinity Q(Ai,0) 7^ Q{Ai — Ai,0) so that the 
single-current fluctuation theorem does not hold in gen- 
eral although the two-current fluctuation theorem always 
does. Furthermore, we notice that the effective affinity 
A\ = 1.6319 is much smaller than the affinity determined 
by the reservoirs: A\ = /3(/J,xl — /Urn) = 25. 

In Fig. IU the single-current generating function is de- 
picted for the smaller value of the Coulomb repulsion 
/3U = 10 and T 2s /Ti s = 2. Here, the effective affinity 
takes a larger value, but again the asymmetry of the gen- 
erating function is still manifest. The shape of the gen- 
erating function now deviates from the parabolic shape 
seen in Fig. [3]as its maximum approaches the unity value. 

Although the ratio of the rate constants is of order 
unity in both Figs. [3] and 0] the difference between the 
generating function and its symmetric is smaller than 5% 
and could remain unobservable if the counting statistics 
was not precise enough. 

G(M) 



0.05 




FIG. 3: The cumulant generating function versus the counting 
parameter Ai at A2 = and the symmetric function with 
respect to the effective affinity Ai = 1.6319 (dotted-dashed 
line) for the parameter values f3U = 30, /3ei = 0, /3e2 = 35, 
0I1XL = 25, PiJ, m = 0, = 70, /3/i 2R = 0, r 1L = r m = 
Til = Tir = 1, T2L = T2R = T2L = ^2R = 1. 
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FIG. 4: The cumulant generating function versus the counting 
parameter Ai at A2 = and the symmetric function with 
respect to the effective affinity Ai = 16.8356 (dotted-dashed 
line) for the parameter values (3U = 10, /3ei = 10, /3t2 = 35, 
Phxl =_ 25, /3/xm = 0, Pfj, 2h = 70, Pjl 2 r = 0, T 1L = r m = 
Til = Tir = 1, F 2 l = I^r = T2L = r 2 R = 2. 




FIG. 5: The cumulant generating function versus the counting 
parameter Ai at A2 = for different values of the electrostatic 
coupling parameter PU. The other parameters take the values 
Pei = 10, /3e 2 = 35, Pjxi h = 25, Pfim. = 0, /3/i 2L = 70, 
l3fi2R = 0, Til = Tir = Til = Tir = 1, T 2 l = T 2 r = 100. 

Figure [5] shows the deformation of the generating func- 
tion Q(Ai,0) as the electrostatic coupling parameter U 
varies from zero to j3U = 20 for r2 S /ri s = 100. In 
the absence of electrostatic coupling, the single-current 
fluctuation theorem holds in the circuit No. 1 since it is 
decoupled from the rest of the system. In this case, the 
affinity takes the value A\ = 25 determined by the two 
reservoirs of this circuit, as seen in Fig. [5] However, the 
non-trivial root A\ of the generating function decreases 
as the Coulomb repulsion U increases, showing the back- 
action effect of the secondary circuit due to the capacitive 
coupling. In the same progression, the maximum of the 
generating function is also reduced. 

For the ratio of rate constants taken in Fig. [SJ the 
generating function is already practically indistinguish- 



able from its symmetric Q(Ai — Ai, 0) so that the single- 
current fluctuation theorem is already effective and the 
considerations of Sec. [V] apply. In particular, the effec- 
tive affinity is now very well approximated by Eq. (|83p . 



D. The large current ratio limit and the effective 
affinity 

In the limit where the ratio of rate constants tends to 
infinity, the generating function becomes identical with 
its symmetric, as argued in Sec. |V] In order to verify this 
prediction, we depict in Fig. [6] the difference between 
both functions versus the counting parameter Ai. We 
observe in this figure that the difference is reduced by one 
order of magnitude each time the ratio of rate constants 
T2/T1 is increased by the same factor. Consequently, 
the single-current fluctuation theorem is well established 
in the large ratio limit T2/T1 — > 00. In this limit, the 
effective affinity is given by Eq. (l83l) . 

lnteOj.oj-edj-vo)] 




FIG. 6: The difference between the cumulant generating func- 
tion and its symmetric with respect to the effective affinity Ai 
versus the counting parameter Ai at A2 = for the parameter 
values PU = 10, /3ei = 10, /3e 2 = 35, P^ = 25, Pum = 0, 
/3^2L = 70, PH2R = 0, Til = Tm = r 1L = r m = 1, and T 2 = 
r 2L = T 2 r = r 2L = r 2R = 1, 10, lOO, lOOO, lOOOO. As ob- 
served in Figs. Oandgl the difference Q(Ai, 0)— Q(Ai — Ai, 0) 
is positive for Ai < Ai/2 and negative for Ai > Ai/2. Here, 
we only depict the difference for Ai < Ax/2. The other half 
has a similar structure if the absolute value of the difference 
is taken before the logarithm. 

The effective affinity is depicted in Fig. [7] as a function 
of the energy pa of the QD No. 1 for 0U = 15. We 
observe that the effective affinity takes the actual value 
([51)) determined by the reservoirs for either low or large 
values of the energy /3ei. However, the effective affinity 
undergoes a significant reduction in between, down to a 
minimum of about A\ ~ 0.45 x A\. The function has a 
characteristic shape, which can be explained in terms of 
the Fermi-Dirac distributions entering in the expression 
of the effective affinity. Supposing that /12R < 62 < 
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fJ-2L~U and < U < Mil — MiRj we nn d that the effective 
affinity is approximately given by 



A, 



/3(Mil - Mm) 
Pi-et - U + hxl) 
/3(mil - MiR - U) 
P(ei ~ Mm) 
/3(mil - Mir) 



for ei < Mir — C7 
for /L«i R - J7 < ei < /iiR 
for /L«i R < ei < /L«il - U 
for /LtiL - U < ei < /LtiL 
for ^il < ei 

(94) 

up to corrections that are smaller than /3 = (k^T) 1 
in the zero temperature limit T — > 0. Crossovers hap- 
pen where the energy e\ coincides with the values of the 
chemical potentials of the left- and right-hand reservoirs 
and the chemical potentials reduced by the Coulomb re- 
pulsion U. The slope of the effective affinity versus Pt\ 
is successively {0, — 1, 0, +1, 0}, as seen in Fig. [71 Ac- 
cording to Eq. (|94p , the minimum value of the effective 
affinity is approximately given by A\ ~ A\ — f3U — 10 in 
the middle interval /3/iiR = < /3e\ < /3mil — PU = 10. 
The affinity A\ = 25 of the reservoirs is recovered for 
/3ei < /3mir ~ PU = -15 and for f3ei > /3mil = 25, 
which explains the features observed in Fig. [7] 





FIG. 8: The effective affinity J83JI of the QD versus the di- 
mensionless energy fie\ of its level for the parameter values 



is approximately given by 



- Mir) 
/3(-ei - C/ + mil) 


/3(ei - Mir) 
/3(mil - Mir) 



for ei < Mir - U 
for /iiR - [7 < ei < /iiL - C7 
for mil - ^ < ei < Mm 
for MiR < £ i < Mil 
for /xil < ei 

(95) 

up to corrections that are smaller than j5 — (k^T) 1 in 
the zero temperature limit T — > 0. In the middle interval 
/?Mil — PU = —7.8 < /3ei < /3mir = 0j the minimum 
effective affinity reaches a value that vanishes in the low 
temperature limit T — > 0. The actual value of the affinity 
Ai = 25 is recovered for /3ei < /?mir — PU = —32.8 or 
pe x > /3/ii L - 25. 



-40 -20 20 40 r 

FIG. 7: The effective affinity J83J of the QD versus the di- 
mensionless energy j3e\ of its level for the parameter values 
/3U = 15 and P? |) -(P )) , 

Equation (|94p predicts that the minimum value of the 
effective affinity could be further decreased by increasing 
the Coulomb repulsion U. This is indeed the case as ob- 
served in Fig. [8j which depicts the effective affinity versus 
the energy t\ now for the value (|56| . Here, we see that 
the effective affinity may vary from the maximum value 
given by the affinity A\ = 25 imposed by the reservoirs 
down to the very small minimum value A\ ~ 0.083565 
at f3ei ~ —3.9525. In particular, the value A\ ~ 3 which 
has been experimentally observed in Ref. [H] is reached 
for 0ei ~ 2.2. 

If the condition /i2R < £2 < M2L — U is still satisfied for 
the parameter values of Fig. [5J the Coulomb repulsion 
is now larger than the difference of chemical potentials: 
U > Mil — MiR- hi this other regime, the effective affinity 




FIG. 9: The effective affinity (f83|) of the QD versus the dimen- 
sionless electrostatic coupling constant f3U for the parameter 
values Pei = and pT fl -lp )! . 

The dependence of the effective affinity (1531 on the 



12 



Coulomb repulsion is shown in Fig. [9] for a given value 
of the energy j3e\ — 0. Here also, the effective affinity 
can be reduced down to a much lower value than the one 
determined by the reservoirs. By a reasoning similar to 
the one used to get Eqs. (|94|) and ([95]) . we can obtain 
the approximate dependence of the effective affinity on 
the parameter U under the conditions /^ir — t\ < < 
Mil - ei < M2L - £2 as follows: 

' 0(-U + Mil - Mm) 

for < U < Mil - £i 
/3(ei - Mm) 
^ I for mil - ei < U < M2L - £2 

1 ~ I P(U + £ 2 - M2L + ei - Mir) 

for M2L - £2 < U < M2L ~ £2 + Mil - £i 
/3(Mil - Mir) 

for fi 2 L - £2 + Mil -e x <U 

(96) 

up to corrections smaller than /3 = (k-gT) 1 as T — > 0. 
The piecewise linear approximation obtained from the 
Fermi-Dirac distributions here also explains the succes- 
sive slopes —1, 0, +1, and 0, observed in the plot of the 
effective affinity versus j3U. We notice that the differ- 
ent linear pieces of the approximation match together at 
the crossover values of the variable /3U. The minimum 
value is reached in the interval /3(mil — e i) = 25 < j3U < 
/3(M2L — £2) = 35 while the affinity Ai = 25 of the reser- 
voirs is recovered for f3U > /3(M2L — £2 + Mil — £1) = 60, 
as indeed confirmed by Fig. |H1 

The lowering of the effective affinity under specific con- 
ditions can be explained in the present model as the ef- 
fect of the back action of the secondary circuit interact- 
ing with the observed quantum dot. The charging and 
discharging rates of the quantum dot can be drastically 
modified by the coupling to the secondary circuit. In 
this way, the effective affinity can be much reduced in 
some regimes which are determined by the value of the 
energy e\ of the quantum dot with respect to the values 
of the chemical potentials and the electrostatic coupling 
parameter U . This back-action effect tends to disappear 
as the temperature increases at constant voltages. 

We emphasize that the reduction of the effective affin- 
ity is not caused by the Coulomb drag. Indeed, there is 
no Coulomb drag for the conditions chosen in the present 
section because we have here taken rate constants such 
that Tj = Tj as in Eqs. fl55J} and ©. Therefore, the 
Onsager coefficient (IT7)) vanishes together with higher- 
order coefficients according to Eqs. (|48|) and (|49|) and 
the Coulomb drag does not manifest itself for the condi- 
tions taken in Figs. 



VII. CONCLUSIONS 

In the present paper, we have reported the study 
of the single-current fluctuation theorem in a Hamilto- 
nian model of quantum electron transport in two capac- 
itively coupled channels, each containing a quantum dot 



(QD). 22 Such a system is similar to the electronic devices 
used in typical counting statistics experiment s 18 ' 31 where 
the current in one circuit can continuously monitor the 
state of the QD in the other circuit thanks to the capac- 
itive coupling. The model allows us to investigate the 
effects of the back action of the monitoring circuit on the 
counting statistics in the light of the so-called fluctuation 
theorems. 

Since both circuits are capacitively coupled and mi- 
croreversibility holds for the total Hamiltonian , a fluc- 
tuation theorem is satisfied for the two currents flowing 
across the system. This two-current fluctuation theo- 
rem (j59"| or pTj) relates the counting statistics of op- 
posite random electron transfers in both circuits to the 
affinities or thermodynamic forces (|55")) - (l56l driving the 
system away from equilibrium. The fluctuation theorem 
is valid far from equilibrium in the strongly nonlinear 
regimes encountered in electronic circuits composed of 
quantum dots and quantum point contacts. 

However, in counting statistics experiments, one cir- 
cuit is used to monitor the current fluctuations in the 
other circuit so that the counting statistics cannot be car- 
ried out on both currents together and is thus restricted 
to a single current. Accordingly, such experiments can 
only test a single-current fluctuation theorem. In general, 
the two-current fluctuation theorem does not imply the 
single-current fluctuation theorem except under certain 
conditions 2 ^ or in some limits as we have demonstrated 
in the present paper. 

In Sec. HVi we have studied the limit of large capac- 
itive coupling between both circuits. In this limit, the 
state of simultaneous occupancy of both QDs in the two 
parallel channels is at a so high energy that it is energet- 
ically forbidden. The consequence is that the two single- 
occupancy states are separately accessible only from the 
empty state and the single-current fluctuation theorem 
holds with respect to the affinity determined by the elec- 
trochemical potentials of the reservoirs. 

In Sec. El we have instead considered the limit where 
the current in one circuit is much larger than in the other 
circuit. Indeed, a large current ratio is a key feature 
of typical counting statistics experiment s 18 ' 31 where the 
current ratio reaches values as high as 10 7 -10 8 . The cir- 
cuit with the very large current performs the continuous- 
time monitoring of the quantum state of the QD in the 
other circuit. In such a limit, the charging and discharg- 
ing rates of the slow QD take values averaged over the 
very fast fluctuations of the monitoring circuit. This 
is the essence of the back action of the monitoring cir- 
cuit onto the QD circuit. As a consequence of the large 
current ratio limit, the single-current fluctuation theo- 
rem holds but with respect to the effective affinity (|83[) . 
which can significantly differ from the actual value of the 
affinity determined by the reservoirs of the correspond- 
ing circuit. This modification of the affinity is due to 
the capacitive coupling between both circuits, as shown 
in particular by Eq. (|94[) . In terms of the parameter 
U of the Coulomb electrostatic interaction appearing in 
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the Hamiltonian ([I]), the affinity is lowered according to 
A\ ~ Ai — j3U under the conditions specified around Eq. 
(|94|) . This result explicitly expresses the effect of the 
back action between both circuits on the single-current 
fluctuation theorem. This back-action effect can be re- 
duced if the Coulomb repulsion U is decreased, but the 
monitoring circuit can no longer resolve the two states 
of the QD as in Fig. [5] if U is too small. On the other 
hand, the back-action effect is also reduced for large val- 
ues of the Coulomb repulsion as shown by Eq. (I96|) and in 
Fig. [S] Indeed, for a large Coulomb repulsion, the affin- 
ity recovers the value determined by the reservoirs and 
the back-action effect disappears. This case corresponds 
to the situation considered in Ref. [34] where a quantum 
fluctuation theorem has been obtained in a multiple mea- 
surements scheme. 

From a general viewpoint, the two-current fluctuation 
theorem implies the non-negativity of the entropy pro- 
duction in agreement with the second law of thermody- 
namics. The dissipation of energy can thus be evalu- 
ated in the electron transport process used to perform 
quantum measurement in the experiments of Refs. [T8ll3ll . 
This dissipation of energy accompanying quantum mea- 
surement is expected on fundamental ground.— The ne- 
cessity of resolving the QD state in real time has for 
direct consequence that the dissipation in the monitor- 
ing circuit is much higher than in the QD by a factor 
Hz/TIi = (A2/A1) x (J2/J1) of the same order of mag- 
nitude as the current ratio J%/Jx. If the QD state is 
monitored with a sampling time At, the secondary cir- 
cuit playing the role of the detector should have tran- 
sitions on equal or shorter time scales according to Eq. 
(|63|) . Since the secondary circuit is driven out of equi- 
librium by the affinity A2, its electron current should 
satisfy J 2 > (Ai) -1 , so that the dissipated power should 
be bounded by U 2 = k B TA 2 J 2 > k B TA 2 (At)- 1 . The 
higher the time resolution, the higher the dissipation 
rate. 

In summary, we have shown that the single-current 
fluctuation theorem is valid under different limiting con- 
ditions and provided a fundamental understanding of the 
back-action effect of the monitoring circuit on the affin- 
ity of the monitored circuit, as observed in Ref. 0. The 
present study extends the analysis of Ref. IT9II21I in show- 
ing how the effective affinity of the single-current fluctu- 
ation theorem can be directly expressed in terms of the 
parameters entering the Hamiltonian of the system. Sev- 
eral issues are left open such as the facts that a quantum 
point contact has specific transport properties and that 
the counting statistics is perf ormed by two QDs in series 
in the experiments of Ref.ll8l We hope to report on these 
issues in a forthcoming publication. 
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Appendix A: Calculation of the average currents 

In this appendix, two different methods are developed 
in order to calculate the average currents given by Eq. 

in terms of the leading eigenvalue Q of the eigenvalue 
problem (|30|) of the four- by- four matrix (|31l) . 

The first method starts from the eigenvalue equation 
(|30]) for the right eigenvector v associated with the eigen- 
value Q and from the adjoint equation 

L T -u=-Qu (Al) 

for the left eigenvector u where T denotes the transpose 
of the matrix. The left and right eigenvectors satisfy the 
normalization condition 

u T • v = 1 (A2) 

Accordingly, the eigenvalues is given by 

Q = -u T L v (A3) 

Taking the partial derivative d a with respect to the 
counting parameter A Q of Eqs. (|A2P - (|A3j) and using Eq. 
(IBTfl) and Eq. (IA1[) . we obtain the following expression 
for the average current: 

J a = d a Q\ x=Q = -u T • d a L ■ v| A=Q (A4) 

Since the left and right eigenvectors are given at A = 
by 




in terms of the probabilities ([27)1 . we get Eqs. (j43")) and 
(|44l) for the average currents. 

With the second method, the average currents as well 
as the linear response coefficients are directly calculated 
in terms of the characteristic determinant (|3"7j) of the 
four- by- four matrix (|31[) . This determinant is a poly- 
nomial of fourth degree: 

Q 4 + C3Q 3 + C 2 Q 2 + C X Q + Co = (A6) 

where the coefficients depend on the parameters of the 
model as well as on the counting parameters Ai and A2. 
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We notice that the last coefficient is just the determinant 
of the matrix (f3"Tj): C = det L. 

Since the matrix (|31j) reduces to the matrix of a jump 
stochastic process conserving probability if Ai = A2 = 0, 
the leading eigenvalue vanishes in this limit: 



0(0,0) = 



(A7) 



Since the average currents are given by Eq. (|42|) , we take 
the partial derivative d a of the characteristic determinant 
with respect to the counting parameter A Q to get 

(4Q 3 + 3C7 3 Q 2 + 2C 2 Q + Ci) d a Q 

+ d a C 3 Q 3 + d a C 2 Q 2 + d a Ci Q + d a C Q = 

(A8) 

Now, the counting parameters must be set equal to zero 
and, according to Eq. (|A7|) . the average current is thus 
given by 



Jot 



d a Co 



(A9) 



Using the symbolic manipulation software 
Mathematical we can evaluate the derivative d a Co if 
the reservoirs of the circuit a are at equilibrium, i.e., 
if its electrochemical potentials are equal so that its 
affinity is vanishing: A a — /3(/x q l ~ Mor) — 0. The 
result is that this quantity vanishes under the condition 
Tj = Tj even if the other circuit is out of equilibrium, 
which establishes Eq. f|48|) . 

We notice that the Onsager coefficient can also be ob- 
tained in the same way. Using Eq. (|4^|) and taking a 
further derivative of Eq. (IA8|) with respect to the other 
counting parameter A^, we get 



L 



d a dpC 



2Ci 



(A10) 



A=0,A=0 



which is used to obtain Eq. (j47|) with the symbolic ma- 
nipulation software Mathematical 



A=0 
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